Vortex pinning by meandering line defects in planar superconductors 
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To better understand vortex pinning in thin superconducting slabs, we study the interaction of 
a single fluctuating vortex filament with a curved line defect in (1+1) dimensions. This problem 
is also relevant to the interaction of scratches with wandering step edges in vicinal surfaces. The 
equilibrium probability density for a fluctuating line attracted to a particular fixed defect trajectory 
is derived analytically by mapping the problem to a straight line defect in the presence of a space 
and time-varying external tilt field. The consequences of both rapid and slow changes in the frozen 
defect trajectory, as well as finite size effects are discussed. A sudden change in the defect direction 
leads to a delocalization transition, accompanied by a divergence in the trapping length, near a 
critical angle. 
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I. INTRODUCTION 

Understanding the physics of thermally wandering 
lines interacting with extended defects is crucial in ex- 
plaining the properties of planar high Tc superconduc- 
tors permeated by columnar pins or twin boundariesi. 
Considerable progress is possible theoretically in planar 
superconductors, when platelet samples are permeated 
by a single sheet of vortex lines due to a small in-plane 
magnetic fielc&^. Related problems arise on vicinal sur- 
faces, where thermally fluctuating step edgea^ can inter- 
act with scratches, grain boundaries or terraces created 
by lithography. In superconductors, a zoo of different 
systems whose properties are determined by the inter- 
play of the interacting vortex filaments, thermal fluctu- 
ations and pinning has been thoroughly explored when 
the attracting defect is straight or vortices follow trajec- 
tories determined by point disorderi2i^*S*iSi. However, 
little work has been done so far when the defects them- 
selves follow a controlled but non-trivial trajectory. On 
vicinal surfaces, such defects could be created by etching 
a wavy, semipermanent scratch, or studied in the con- 
text of the square terraces created lithographically by 
Lee and BlakelyS. One might also consider a thin high 
Tc superconductor sample where a wavy notch has been 
etched, although curved extended defects can arise nat- 
urally in the form of grain boundaries in polycrystalline 
platelet superconductorsi^. Finally, close intersection of 
splayed columnar defect &^O'^" ' provide preferred vortex 
tracks with a sudden change of direction. 

In this work we investigate a single thermally fluctuat- 
ing line that interacts with a single quenched meandering 
defect. This situation mimics the dilute flux line limit of a 
planar superconductor with a single pinning defect when 
the vortices are far apart enough to be considered nonin- 
teracting. Although we use the terminology of flux lines 
throughout, similar results should apply to scratch-step 
interactions on vicinal surfaces. The case of many step 
edges (or vortex lines) interacting with a single curved 
line defect will be treated in a future publication^Si. 

In the spirit of the treatment of high Tc superconduc- 
tors in Refi^, we describe the trajectory of a flux line in 
(1+1) dimensions as a classical elastic string x{t) sub- 




FIG. 1: Thermally excited fluctuating string (vortex line or 
step edge) interacting with a meandering attractive line-like 
potential. The dark line indicates the elastic string trajectory 
x(r), the energy of which is given by Eq.Q. 



ject to thermal fluctuations. We assume high enough 
temperatures and sufficiently clean samples so that point 
disorder can be neglected. In a quantum analogySi^, (see 
below), the spatial direction labelled by r plays the role 
of imaginary time. We assume that overhangs are im- 
probable, so the choice of a single valued function to de- 
scribe a string-like flux line is adequate. The trajectory 
of the defect itself (a notch, etched on the thin super- 
conducting slab, say) will be represented (r), and 
attracts the vortex with a short range attractive potential 
V[x-Xo{t)], (see Fig.Hl. 

For the energy of a flux line induced by a magnetic field 
along r in (1+1) dimensions, bounded by time- like initial 
and final coordinates ti and t/, we take the integral: 



E[x{t)] 



dr 



dx{T) 
dr 



+ V[xiT)-Xo{T)] 



, (1) 
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where 7 is the coarse-grained line tension. The finite 
temperature behavior of the string can be determined 
from the partition function, which is a functional integral 
over all possible vortex configurations x{t), 



'Dx{T)e 



'E[x(t)]/T 



(2) 



where we have set Boltzmann's constant, fc^ = 1. 

The statistical physics of the above system is controlled 
by a competition between the two terms of the line energy 
of Eq.(P). The "kinetic energy" (corresponding to the 
line tension of the string) is minimized when the vortex 
line is straight and parallel to the r direction, whereas the 
"potential energy" is minimized when the string resides 
at the deepest point of V[x[t) ~ a;o(T)], namely exactly 
on the attracting defect. Depending on the strength of 
the potential, the stiffness of the string, the temperature 
and the trajectory Xo(t) of the defect, the elastic string 
will either localize near the trajectory of the defect or 
wander far from it. Since in most experimental realiza- 
tions the attractive potential is relatively short ranged, 
in this work we will take V\x\ in Eq.QJ to be an attrac- 
tive delta function, which should adequately describe the 
long wavelength properties of the system. 

The case of Xo(t) being a zeroth or first order polyno- 
mial of r (i.e., a straight line!) has been studied in detail 
in Ref.^, with emphasis on the delocalization transition 
that takes place for large relative tilts of the pining de- 
fect and the flux line. This work adapts the results and 
conclusions of that paper to the case of a "bent defect" 
with a sudden change in direction, and considers more 
general quenched defect trajectories as well. 

This paper will be structured as follows. In Section UTI 
we review the basic mapping from a classical descrip- 
tion of the thermally fluctuating flux line to a (non- 
Hermitian) quantum mechanical problem, and use this 
mapping to analytically compute the partition function 
(i.e., the propagator for the equivalent quantum prob- 
lem) for the straight defect case. Section lllll explains how 
to use this propagator to derive exact expressions for a 
simple realization of a "bent" defect trajectory and de- 
termine how the defect probability distribution responds 
to the presence of the kink. Section Hvl describes a per- 
turbative solution in the "Born-Oppenheimer" limit of a 
slowly varying Xo (r) . Finally, Appendix g contains some 
calculational details, relevant to Sections ITU and ITTll 



II. PARTITION FUNCTION FOR A FLUX LINE 
INTERACTING WITH A LINEAR DEFECT 

In this section we briefly review the connection be- 
tween the classical path integral of Eq. Q and the "quan- 
tum mechanical" formulation, which follows from a trans- 
fer matrix techniquei^. We then proceed to analytically 
compute the partition function that governs the proper- 
ties of the thermally excited vortex line interacting with 
a straight line defect. 



For clarity we denote the partition function of Eq. ^ 
as Tj; Xi, r^; Xo(t)], where Xi and Xf are fixed ini- 

tial and final positions of the flux line. (We may eventu- 
ally choose to integrate over these endpoints.) Since the 
potential is only implicitly dependent on r through the 
defect coordinate Xo(j), we can map the problem by a 
change of variables a;(r) y(j) = x(t) — Xo{t) onto one 
of a straight defect in a r-dependent external tilt fleld 
/i(r) — dxo{T)/dT that imposes tendency to drift on the 
flux line. The Jacobian of the transformation is unity, so 
the measure of the path integral remains the same. After 
this transformation, the partition function reads: 

Zt[xf - Xo{Tf),Tf;Xi - Xo{n),n;h(T)] = J I?y(r)e~"^'^'^"'''^ 

(3) 

where 



E[y{r)] ^JdT 



7 f dyir) 
dr 



+ hiT)\ -Vo5{y) 



(4) 



The sign of Vo > is chosen so that the potential is 
attractive. By construction, the old and new partition 
functions must agree: 

Z[Xf, Tf, Xi, Ti]Xo{T)] ^ Zt[Xf-Xo{Tf), Tf, Xi-Xo{Ti), Ti] /i(t)]. 

(5) 

One way around the technical difficulties (see, e.g., 
Refii^) of the evaluation of the partition function of 
Eq. (PJ due to the singular nature of the delta function, is 
to recast the problem in quantum language. The classi- 
cal partition function then becomes a matrix element of 
the imaginary time evolution operator: 

Zt[xf,Tf-x,,n-h{r)] = (a;/|r,{e'^7''-«M/i^}|x,), 

(6) 

where the time ordering operator Tt is required by the 
time-dependent Hamiltonian: 



n{T) 



'2^ dx^ 



Hr)T 



d_ 

dx 



Vo 5[x). (7) 



The mapping onto imaginary time quantum mechanics is 
clear if we let T — > ?i and 7 — > m, where m is a particle 
mass. The term /i(t) T makes the Hamiltonian non- 
Hermitian. 

Before considering curved defect trajectories, we ana- 
lytically compute the matrix element Zt[x,T;y,0;h] for 
constant tilt h, by expanding with respect to the com- 
plete set of right and left energy eigenstates {\iPk)r} and 
{ L {i/JK I } of the non-Hcrmitian Hamiltonian of Eq. Q : 

Zt[a:,r;z/,0;/i] =^(x|VK)Ke-^-^/^i(V'if|2/), (8) 



K 

with Ek being the K^^ eigenenergy. Since both the 
Hamiltonian and the propagator are r-translationally in- 
variant, without loss of generality we have chosen zero 
as the origin of imaginary time. For notational simplic- 
ity we have also rescaled the spatial coordinates {x,t) 
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and tilt h so that [x] = T^/-/Vo, [r] = 2T^/jV^, [h] = 
Vo/T. After these rescahngs, the evolution equation for 
Zt[x,T]Xi,Ti;h] = Z[x,t\ reads: 



dZ\ 



dr 



-HZ[x, t] 



ydx"^ 



2h{T)—+2S{x)]Z[x,T]. (9) 



As explained in Ref.^, there is a critical value he of the 
tilt, above which the flux line delocalizes from the defect. 
Our rescaled units are chosen so that this value is he = 
1 for our delta function potential. Notice that in the 
rescaled units the track of the straight line defect is given 



dXo{r) 



2h. 



by ST 

The set of normalizable eigenstates of Eq. {T)) is differ- 
ent above and below the critical tilt he (see Appendix^; 
however, the propagator when the lateral dimension 
of the sample tends to infinity has the same form for both 
regimes, namely: 



r 




X 



FIG. 2: Flux line near a kinked defect: The defect, initially 
parallel to the average direction of the freely fluctuating vor- 
tex suddenly tilts at r = n. 



Zt[x,T;y,0;h]^e-''^--y^-'' ^ 
g-lal-kl+T 



2V^ 



erfc ( ^ 



(10) 



where erfc(a;) = f dt e * . The first term inside the 

square brackets is independent of the defect, transla- 
tionally invariant in the x direction and describes the 
random thermal "diffusion with drift" of the flux line, 
Z ~ exp[— (x -I- 2/it)^/4t]^^ . This term dominates for 
small times and/or large |a;|, \y\, i.e. far from the defect. 
The second term becomes important for sufficiently large 
imaginary times or small \x\, \y\. The delocalization tran- 
sition manifests itself in the exponential growth of the 
second term for large \x\, \y\ when the tilt h in the pref- 
actor exceeds the critical tilt he = 1. Upon returning to 
the original units of x and r, we see that the effect of the 
defect becomes apparent for r > r * ~ T^/jV^, provided 
I a; I and \y \ are within a localization length a;* — T'^ /jVo of 
the defect at the origin. For Vo 0, the first term in the 
brackets scales as 1/Vo and will dominate over the second 
term which becomes independent of Vo- This limit leads 
to the expected free flux line propagator. These results 
can be easily translated in the original frame of reference 
(where the defect is tilted) by the use of Eq. jSJ . 



III. PROBABILITY DISTRIBUTION OF THE 
FLUCTUATING VORTEX FILAMENT 

In this section we use the above propagator to ana- 
lytically compute the probability distribution function of 
a flux line interacting with a r-dependent defect. As a 
simple, yet revealing example we study a line defect con- 
sisting of two straight lines joined at an angle, at r = ti. 



namely: Xo{t) = 2h {t ~ ti)9(t — ti), see Fig. |21 Here, 
Q{x) is the step function, Q{x) — 0, x < 0, Q{x) = 1, 
X > 0. Any sufficiently well behaved defect trajectory 
can be approximated by a series of finite segments; the 
method described below can be straightforwardly gener- 
alized to more complicated piece-wise linear configura- 
tions. 

When h < 1, and in particular when h = for the first 
segment of the defect, the system relaxes to the localized 
ground state exponentially fast as t becomes large. A 
ground state initial condition \tp^) = \tpg) at r = is 
thus the experimentally relevant case when the kink at 
Ti is sufficiently far from the boundaries. For t — 
we choose free boundary conditions embodied in the fi- 
nal state('0^| = J dy{y\, meaning that no restriction is 
imposed on the exit point of the vortex from the sample. 
In what follows, unless otherwise stated, we work in the 
thermodynamic limit in the x direction, so integrations 
such as J dx may be understood as dx. 

The probability distribution for the flux line on some 
internal time slice t is formally equal to: 

P{x,t) = l(ijf\S{Lr,r)\x){x\S{T,OW), (11) 



where 



(12) 



is the time ordered evolution operator associated with 
the Hamiltonian 



nir) 



' 27 dx^ 



Vo i{x - Xo{t)\. 



(13) 



In Ea. (|ll|l . Z is used to represent the partition function 
for the entire length of the system on both sides of the 
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kink, 

Z={ijf\S{Lr,QM)- (14) 

Since the defect trajectory Xo{t) is piecewise straight, 
we can spUt the factors in the numerator of Eq. into 



The last propagator inside the integrand effectively lo- 
calizes the vortex near the defect at Xo{t) = case. 
Since the initial condition is the ground state of the time 
independent Hamiltonian, this last factor is just equal 
to e~l^le"^«'^ = e~l^le'^, since the ground state energy 
in our units is Eg ~ —1. The middle propagator in 
Eq.lHsll follows from Eq.lUnil by setting h ^ 0. The left- 
most factor in the integrand can also be expressed in 
terms of Ea. H10(l . if we first perform the tranformation 
of Eq.(|SJ). Then, since y' is a dummy integration vari- 
able and Xo{ti) = 0, it is easy to see that Ea. (|15|l can be 
rewritten in terms of the transformed partition function 
appearing in Eq.© as: 

P{x,t) = ^ ^d2/d?/'Zt[y,i^;y',ri;/i]- 

•Zt[y',ri;x,T;0]e"l-le^ (r < n). (16) 

Although the y-integration can be computed analytically 
(see Appendix^, we have resorted to numerical integra- 
tion to handle the integration over y' . 

Upon following a similar procedure for r > ri, we ob- 
tain: 



P{x,t) = 




dyZt[y,Lr;x - Xo{T),T;h 



■ J dy'Zt[x-Xo{T),T;y',n;h]e-^y'^e^\ 

(r>ri). (17) 

As shown in Appendix „ in this case both the integrals 
can be done analytically. 

Having now computed P{x, t) for all imaginary times, 
we can examine the flux line probability distribution for 
this simple realization of an attractive line defect that 
suddenly changes direction. We have studied two cases, 
one where the tilt of the upper segment of the defect is 
below the critical value he = 1 for delocalization (Fig.O 
and one where it is above (Fig.^. 

Note the following: For time-like coordinates below 
the kink at ti, the probability density does not differ 



two pieces such that dxo/dr = const. For computational 
purposes, the cases r > ri and r < ri need to be treated 
separately. 

First, for r < ti we expand the probability distribution 
in position eigenstates as: 



(15) 



h=0.8 
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FIG. 3: Contour plot for the probability density distribution 
of the line defect for a scratch everywhere below the critical 
tilt angle: pinned regime. The darker areas correspond to 
higher probability. The tilt of the upper segment of the defect 
is /i = 0.8 < he. = 1, the cusp is at n = 7/3 ~ 2.33 and the 
total length of the sample is hr = 7. The inset shows the 
probability density distribution for: r = 1 solid line, r = 2.33 
dashed line, r = 3 dash-dotted line, r = 4 dotted line. Note 
that the distribution, while remaining localized, develops an 
asymmetry for r > 7/3. 



much from the ground state probability density of Ref^ 
P{x, t) = e~^l^l. However, as the position ti of the kink 
in the scratch is approached from below, there is a broad- 
ening of the probability density on the right or "concave" 
side of the defect, and a similar narrowing on the left or 
"convex" (see Fig. This phenomenon is not directly 
related to the delocalization of the flux line that occurs 
for tilts greater than the critical tilt, because it is present 
for both h > 1 and h < \ (recall that the critical tilt is 
he — 1). For h <1, after this small shifting of the prob- 
ability density in the vicinity of the kink, the probabil- 
ity again resumes the ground state distribution. When 
r ^ Lt- we notice a small broadening of the probability 
which is a boundary condition effect. 



J 



P{x,t) = 1 j j Aydy'{y\S{Lr,T,)\y'){y'\S{T,,T)\x){x\S{TM4'g) 

= ^ JJ '^y'^y'^ly^Lr^y' ,Ti;XoiT)]Z[y' ,Ti;x,T;0]{x\S{T,0)\i)g) 

I 
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h=1.5 




-4 -2 2 4 6 8 10 

X 



FIG. 4: Contour plot for the probability density distribution 
of the defect line whose upper part is above the critical tilt 
angle for the scratch: depinned regime. The tilt of the up- 
per segment of the defect now is/i=1.5>/ic = l, while 
Ti — 2 and Lt — 6. The inset shows the probability density 
distribution for: r = 1 solid line, r = 2.33 dashed line, r = 3 
dash-dotted line, t — 4 dotted line. Note that the probability 
distribution begins to delocalize for r > ri = 2. 



When the tih of the second segment of the defect is 
above the critical tilt he = I (see Fig. 0J), the proba- 
bility density remains localized around the defect only 
for some small distance above the kink. Two local max- 
ima of probability density appear, one localized on the 
defect and the other representing an approximately gaus- 
sian probability distribution. As we move further away 
from the kink, the probability of observing the flux line 
close to the defect reduces, and the weight associated 
with the delocalizing diffusive part increases. 

To better understand the depinning transition as the 
slope of the kink angle approaches the critical value he 
from below, we have computed the average position of the 
fluctuating line {x)r = J dx x P{x,t) for several tilts and 
sample lengths. In Fig. [S] we have plotted the average 
position of the fluctuating vortex line for several total 
lengths in the imaginary time direction Lr- The position 
of the kink is held fixed at ri = 2 and the tilt is h = 0.8. 
We observe that as Lr grows, the distance m of the exit 
point of the vortex line {x)l^ from the defect, indicated 
in the figure, tends to a constant. We have computed 

m{h)= \im {{x)l^ -Xo{Lr)), h < 1, (18) 

by obtaining values for up to L^- — 50 and then extrapo- 
lating Lr to infinity. The convergence becomes very slow 
as h — > 1^. The limiting offset m{h) serves as an indi- 
cator of the depinning transition for h h~ . Similarly, 
for /i > 1, as can be seen from Fig|Sl a corresponding 
indicator is the offset of (x) from the origin, i.e. the 



h=0.8 




-2 2 4 6 81012 



FIG. 5: Average position of the vortex {x)r : pinned regime. 
Here ri = 2 and h — 0.8. The dashed line is is the scratch 
trajectory and the solid lines represent the thermal average 
{x)t for different sample lengths. The arrows mark the aver- 
age distance of the step from the exit point of the scratch. As 
Lr increases, this distance approaches a constant shift in the 
trajectory. 

h=1.4 




0.5 1 1.5 2 

<x>, 

FIG. 6: Average position of the vortex {x)t- depinned regime. 
Here ri=2 and h=1.4. The arrows represent the deflection of 
the vortex trajectory due to the tilted scratch and the dashed 
line is again the defect trajectory, (x) converges to ~ 2.2 
as Lr increases. 

original location of the defect. We now measure a shift 
m{h)= \im {{x)lt -Xo{Q)), h>l, (19) 

which contributes to the total magnetization, in the vor- 
tex line case. By a fitting of m{h) and m{h) to l/\h — 
(see Fig. U)) we find the critical exponent of the delocal- 
ization transition to be i/ = 1.0^^. An identical exponent 
describes the delocalization transition studied in Refi^. 
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>< 30- 




Each of the summations embodied in ^ runs over 

l/'l,l/'2,---,l/'N 

the entire spectrum, while |V'o(0))i? = |V's(0))i?- 

In the hmit e — > 0, we can expand |V'fc-i(ifc — e))_R as 
\il^k-i{tk - e))ii ~ \'4)k-i{tk)) R - €^\ipk-i{tk))R, so: 



L(V'fc(ife)|V'fe-i(ife-i))fl = '^•0fc,v>fc-i-e^Vfe.V'fe-i(*fe)+O(e^), 

(23) 

where 



= {■>Pk{tk)\j^\A-i{tk-i))R.. (24) 



To Hnear order in e, we can rewrite (|23f) as a sum of two 
mutually excluding terms: 



FIG. 7: Shift in the line trajectory near the depinning transi- 
tion; curve to the left of hc=l represents m{h) while that to 
the right represents m{h). 



(25) 



IV. SOLUTION FOR A RANDOM 
TRAJECTORY IN THE ADIABATIC LIMIT 

In this section we use the time dependent perturba- 
tion theory developped in Mostafazadebi^ to analytically 
study the case of a slowly bending defect, i.e. | ^ | ^ 1- 
We will compute the probability distribution defined by 
Ea. (|ll|l . assuming that the initial and final states are the 
ground state, which should be a good approximation for 
physical systems with h =const for r everywhere outside 
the interval (0, Lr). As in Sec. In] we will preform a defect 
coordinate tranformation 



X X + Xo(t) 



(20) 



and reduce the problem to a straight defect in a time 
dependent external tilt field. 

We will assume that |/i(t)| < 1. For any fixed /i(r) we 
can then be sure of the existence of a bound state. We 
suppose for the moment that < oo, which allows us to 
consider a discrete spectrum of free states. We can then 
write the matrix element of the time evolution operator 
as a product over time slices: 



(x|5(r,0)|^3)= hm ^ 



- 5 E^.{t,)e 



X Y[{x\lbNiT))RL{i^k{tkMk-litk-l))R (21) 



N 



fe=l 



where e = t/N, tk = k e and \ipk{tk))R and L(V'fc(ifc)| are 
respectively the instantaneous right and left eigenstates 
of 7i(T) for h = h{T = t^), namely: 



Upon inserting Eq. 1)25(1 back to the matrix element 1(21(1 
and keeping terms only up to the first order in jA^^^^^j 
(which is small if /i(t) is slowly varying), we have: 



'A'^g{T))R - ^ {x\i}{T))R- 

dte-.A^'^*'^^(*')A^,V'.W)- (26) 



n{T)\Mtk))R = E^„{tk)\Mtk)) 



R- 



(22) 



Similarly, the other matrix clement in Eg. ((11(1 . 

L{iljg{Lr)\S{Lr,T)\x), reads: 

i(^<,(L,)|5(L„T)|.T) = e-.''^ 

'i(^g(r)|a;) - ^ i(^(T)|a 

^'^'dte-/^*'^*'^^(*')A^^,^(i)), (27) 

where AE{t) = E^{t) ~ E^^{t) -f A^^^{t) - A^^^^^{t). 
Upon sending to infinity, the sum becomes an in- 
tegral over wavevectors k and {x\^-)k{T))R becomes the 
eigenstates tabulated in Ea. l(A.2() of the Appendix, with 

h = h{T). 

Provided ^.(t) < 1, to zeroth order the probability 
density is independent of h, P{x,t) ~ e"^'^'. The low- 
est order adiabatic correction to the probability density 
AP{x,t) = F(x,t) - Po(t) reads: 
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AP{x,t) = e 



dk 



27r(l + c) 



- |x| 4'/i(T)a; 



dk 



27r(l + c) 




FIG. 8: Numerical solution AP{x, t) for a slowly moving de- 
fect. The thick solid line in the center represents the trajec- 
tory of the defect, Xoij) = 18(1 — cos(tjr)), where ui = 0.1. 
The time domain shown corresponds to slightly less than 5% 
of a complete period. The dashed and solid lines are respec- 
tively negative and positive contours of the adiabatic correc- 
tion AP(x,r). On the right: snapshots for times: r = 0.03 
and r = 1.50 of A.P(x,t) for both the analytic and the nu- 
merical solution. 



FIG. 9: Numerical solution A.P{x,t) for a rapidly moving de- 
fect. The thick solid line in the center represents a half-period 
of the defect trajectory, Xo(r) = 3.6(1 — cos(a;T)), where 
LJ — 0.5. The dashed and solid lines are respectively negative 
and positive contours of the adiabatic correction AP{x,t). 
On the right: snapshots for times: 1.02, 4.02 and 5.88 of our 
approximate solution of AP{x, r) compared with the numer- 
ical result. 



with c = — (see Appendix,). J AP{x,T)dx — 

and P{x, t) is normalized to unity. 

Upon substituting x ^ x — Xo{t) we can trivially tran- 
form back to a "moving defect" frame of reference. In this 
frame, the zeroth order term reads Pair) = e~^l^~^°('^)l . 

To check the accuracy of this result, we have used 
the moving defect frame and numerically integrated the 
imaginary time Schroedinger equation for a slow sinu- 
soidal defect trajectory Xo{t) — 18[1 — cos(O.lr)] with 
3 to obtain AP{x, r) = P{x, t) ~ e-2b-^o(r)| ^ r^Yic 
results are shown in Fig. |HJ^. On the right hand side of 
FiglSlwe plotted two snapshots of the theoretical (lowest 
order correction) and numerical AP(a;, r) for two differ- 
ent times. The agreement between the two results is 
fairly good, even for x and t close to the boundaries. 
Note that the probability density piles up on the side of 
the radius of curvature, and is depleted on the opposite 
side. 

However, the agreement deteriorates for rapidly chang- 
ing defect trajectories. In Fig. |51 we have plotted the 



numerically evaluated AP{x,t) for Xo{t) ~ 3.6[1 — 
cos(0.5t)] with Lt ~ 6.3 and compared with the ana- 
lytical result. The three snapshots on the right indicate 
significant deviation from the theory in the center of the 
range of r-values, where ^ is changing rapidly. Note 
from the more exact numerical solution that probability 
again piles up on the side of the radius of curvature. 



Acknowledgments 

We would hke to thank M. Aziz, J. Blakely, G. Chan, 
V. Ignatesku and E. Williams for useful conversations. 
This work was supported by the National Science Foun- 
dation through NSF Grant DMR-0231631 and the Har- 
vard Materials Research Laboratory through NSF Grant 
NMR-0213805. 



8 



APPENDIX: ANALYTICAL COMPUTATION OF 
THE PROBABILITY DENSITY 

In this Appendix we derive a number of results given 
in the main part of this paper. In most of what follows 
we assume that the lateral dimension of the planar su- 
perconductor Lx becomes infinite, although we shall be 
interested in the dependence of physical quantities on the 
size of the timelike dimension L^- All the quantities are 
given in the dimcnsionlcss units discussed in SecHll 



Above the delocalization transition, the nature of the 
eigenstates changes and there is no longer a bound state. 
The non-Hermitian nature of the Hamiltonian |2J) now 
becomes important. The right eigenstates for h > 02i 
are (see Ref,^): 



ikx 



ce 



ikx — 2hx 



■ c e 



and the left ones: 



for a; > 
for a; < 0, 



(A.2) 



1. Eigenstates of the Hamiltonian 

The well known ground state of the Hamiltonian (jjj) 
with /i = (in rescaled units) is simply il)g{x) = e"'"^' 



with eigenenergy Eg 



-1. 



20 



The odd and even 



parity extended eigenstates are respectively ^'^'^'^{x) ~ 
^i+k'^ (A: cos kx — sin fc I a; I) and il)™'^'^{x) ~ sinfca: with 
Ek = k"^ and k chosen to satisfy periodic boundary con- 
ditions in the x direction. 

This complete set of eigenstates can be used in the 
expansion ||HJ), the sum being turned into an integral as 
Lx ^ oo, which allows direct computation of the prop- 
agator. Provided \h\ < 1, the case h can be re- 
duced to /i = via the "imaginary gauge transforma- 
tion" e^-^H{h)e-'''' ^ H{0). In this case, Eq.lUHIl can be 
obtained by straightforward evaluation of the integral: 



Z4a;,r;y,0;0] = e"l^l-l^l+^ + - / & e" 

Jo 



{kcoskx — sin/(;|a;|)(/ccosfc?/ — sinfc|?/|) 



l + fc2 



sin fca; sin A;?/ j . (A.l) 



i'Lkix) 



(1 + c)e'''^ for X > 

^rkx ^ ^^-tkx+2hx for x < 



with a common eigenenergy 

Ek^ k'^ + 2ihk, 

and 



ik — h 

The normalization condition is: 



(A.3) 



(A.4) 



(A.5) 
= 1. 



/ dk\tpk)R L{ipk\ 
Now the propagator is evaluated by computation of the 

integral jT^dk e~'^ ^'ijj^k{x)^^k{x)- Remarkably, the 
final result for the propagator is again Eq. l(Tn|l . 



2. Exact probability distribution of the fluctuating 
string 

In this subsection we sketch details for obtaining re- 
sults used in Sec lIIII 

The y- integration in Eq. can be computed analyt- 
ically and is found to be: 



J dyZt[y,Lr;y',Ti;h] = 1 - 



^t-\y'\+hy'-h''t 



Jiy' ( ph\y'\ 



1 



-erfc 



e~h\y'\ 
1 - h 



-erfc 



2^/^ 



(A.6) 



where t = Lr — ti. 

The y integration in Ea. (|17|l is formally the same as 
the integral IjA.GI) . ti and y' having to be substituted 
with T and x — Xo{t) respectively. The y' integration is 
similarly found to be equal to: 

J dy'Zt[X,T-y',Ti-h]e-\y'\ =Q{h,X,T-Ti)+ 

+ Q{-h,~X,r-Ti), (A.7) 



where X = x ~ Xo{t) and 

Q{h,x,t) = ie(i"2'')*'^erfc + (1 - h)Vi^ + 

+ e(i-M^*+(i-'»)klerfc + (1 _ h)Vi) ] . (A.8) 
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Returning to the original, tilted defect frame of reference, 
we see that Z ~ exp[— a;'^/4r] as expected from Eg. (1131 for 

m{h) is related to the trapping length introduced in^ as a 
measure of the offset of the flux line from the average po- 
sition it would assume without the defect (here that would 
be {x)r = because of the boundary condition) and the 
critical exponent computed here agrees with the exponent 
analytically predicted there. 
^'^ The corresponding tilt fleld for this trajectory is ^(r) = 
0.9sin(ajT), since with our rescaled units /i(r) — ^^§^- 
See for example E. Merzbacher, (John Willey & Sons, New 
York, 1998), Chapter 6. 

The h < eigenstates can be derived by considering that 
nHh) = n{-h). 



